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摘要 :， 自 开始 探测 引力 波 以 来 ， 核 心 夫 缩 型 超新星 就 被 认为 是 一 类 重要 的 引力 波 天 体 物 理 源 。 
可 顾 了 广义 相对 论 中 引力 波 的 相关 知识 ， 介 绍 了 核心 拥 缩 型 超新星 在 Fe 核心 拥 缩 阶 段 和 反弹 阶 
Bt, dt Fe 核 旋转 的 角速度 和 对 流 不 稳定 性 对 引力 波 暴 的 影响 ， 并 且 讨 论 了 核心 夫 缩 型 超新星 爆 
发 的 最 终 产 物 一 一 原 中 子 星 的 振荡 理论 以 及 其 产生 的 引力 波 暴 。 最 后 ， 展 望 通过 电磁 波 和 引力 波 
的 协同 观测 ， 将 来 能 够 更 好 地 理解 核心 雪 缩 型 超新星 爆发 的 全 过 程 ， 并 进而 探知 中 子 星 的 内 部 结 
构 。 
X BH dg: 核心 地 缩 型 超新星 ， 引 力 波 暴 ， 中 子 星 ， 振 荡 模 式 
中 图 分 类 号 : P142.6 文献 标识 码 : A 


1 引 8 


引力 波 是 爱 因 斯 坦 ”” 的 广义 相对 论 预言 的 一 种 时 空 涟 满 ， 它 是 时 空 本 身 的 一 种 波动 现 
象 。 根 据 引 力 波 理论 ， 物 质 产生 的 引力 辐射 的 功率 与 物质 质量 四 极 矩 对 时 间 的 三 阶 导 数 的 
平方 成 正比 。 因 此 ， 质 量 四 极 矩 的 变化 越 快 则 产生 的 引力 波 的 强度 越 强 ;关于 引力 波 的 理论 
可 参考 Maggiore 的 著作 ”。 在 地 面 引力 波 探测 实验 中 ， 引 力 波 探 测 器 主要 分 为 共振 棒 引 力 
波 探测 器 和 激光 干涉 仪 引 力 波 探测 器 两 类 吕 ， 其 中 激光 干涉 仪 引力 波 探 测 器 由 于 探测 灵敏 度 
高 、 频 率 范围 大 ， 很 快 便 成 为 主要 的 引力 波 探测 器 。2016 年 2 月 11 日 ， LIGO M Virgo 合作 
组 宣布 ， 他 们 首次 探测 到 了 双 黑 洞 并 合 产生 的 引力 波 信号 ， 即 GW150914 引力 波 事件 ”。 
如 今 ， 随 着 LIGO 和 Virgo 激光 干涉 仪 引 力 波 探测 器 的 升级 改造 ， 越 来 越 多 的 双 黑 洞 并 合 产 
生 的 引力 波 事件 被 探测 到 ””。 因 此 ， 引 力 波 的 发 现 开启 了 引力 波 天 文学 研究 的 新 时 代 。 
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由 于 实验 室 无 法 产生 可 供 探测 的 引力 波 ， 因 此 ， 大 质量 、 高 速 运动 的 天 体 就 成 为 重要 
的 引力 波源 ， 而 且 质 量 越 大 、 天 体 越 致密 ， 产 生 的 引力 波 的 总 能 量 也 越 大 。 如 双 黑 洞 并 合 
事例 ”， 此 外 ， 还 包括 双 中 子 星 绕 转 ””"、 超 新 星 爆 发 ””、 中 子 星 振 荡 ”” 等 现象 。 传 统 
上 ， 实 验 探测 的 引力 波 波源 主要 分 为 致密 双星 的 绕 转 与 并 合 产生 IRT. ERE 
WU. EERS Zu ums sU 7. RAPE T RUPEM BOR E 
fk, BUC YE RAC. Hie AR A. E —Jr I RE JRUR HARE 
中 往往 需要 考虑 到 流体 动力 学 、 中 微 子 输 运 过 程 、 磁 场 等 ， 因 此 很 难 进行 数值 计算 ， 另 一 方 
面 ， 核 心 夫 缩 超新星 爆发 的 最 终 产物 之 一 是 中 子 星 ， 而 对 于 中 子 星 的 内 部 结构 至 今 也 没有 明 
确 的 答案 。 基 于 上 述 两 方面 的 原因 ， 虽 然 核心 声 缩 超新星 爆发 过 程 会 产生 引力 波 信号 ， 引 力 
波 暴 的 精确 预言 仍 具 有 很 大 的 不 确定 性 。 图 中 是 Ott 等 人 模拟 的 旋转 核心 击 缩 产生 的 引 
力 波 暴 波形 示意 图 ， 其 波形 的 细节 特征 仍 有 待 更 多 研究 。 
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注 : 波源 的 距离 是 10 kpc。 


1 引力 波 暴 波形 示意 图 


下 文 我 们 将 着 重 介 绍 核心 替 缩 型 超新星 爆发 过 程 ， 包 括 核 心 替 缩 、 反 弹 以 及 最 终 形成 
中 子 星 阶 段 产生 的 引力 波 暴 ， 并 借 此 展望 电磁 辐射 与 引力 波 辐射 的 协同 观测 的 前 景 和 意义 。 
本 篇 文章 中 所 有 波源 的 距离 均 采 用 D = 10kpc。 
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爆发 时 释放 的 能 量 约 为 10$ J。 质 量 不 小 于 10 Mo 的 前 身 星 ， 其 中 心 部 分 经 过 核 聚变 后 最 终 
成 Fe 元 素 ， 而 Fe 元 素 有 最 大 的 比 结合 能 ， 所 以 不 会 继续 产生 核反应 释放 能 量 以 抵抗 自 引 
力 。 随 着 聚变 反应 的 不 断 进行 ，Fe 核 的 质量 将 不 断 增加 ; 当 Fe 核 的 质量 超过 钱 德 拉 塞 卡 质 
量 极限 时 ， 自 引力 将 超过 电子 简 并 压 ，Fe 核 将 被 进一步 择 缩 形 成 中 子 星 或 黑洞 。 

FEW ARE, 24 Fe 核 中 心 密度 超过 饱和 核 物质 密度 pvuc © 2.8 x 101* g- cm? 时， 将 
会 产生 反弹 激 波 。 早 期 关于 核心 替 缩 型 超新星 的 研究 认为 ， 反 弹 激 波 携带 的 能 量 足 以 使 超 新 
星 爆发 ， 这 种 机 制 被 称 为 瞬时 暴 。 然 而 ， 后 来 研究 发 现 激 波 在 向 外 传播 的 过 程 中 使 得 Fe 原 

子 核 裂 变 ， 激 波 将 消耗 大 量 的 能 量 ， 最 终 不 可 能 导致 超新星 爆发 ”。 超 新 星 爆发 之 后 形成 
的 中 子 星 称 为 原 中 子 星 ， 其 内 部 温度 很 高 且 通 过 辐射 中 微 子 来 使 得 内 部 温度 降低 ， 在 这 过 程 
中 携带 的 能 量 约 为 10%J， 这 当中 约 有 1% 的 能 量 被 原 中 子 星 外 围 物质 吸收 。 研 究 发 现 原 中 
子 星 外 围 物 质 吸收 的 能 量 或 能 成 功 导致 超新星 爆发 ””。 这 个 过 程 中 微 子 向 外 传播 需要 一 
定 的 时 间 才 能 将 能 量 传递 给 外 围 物质 ， 因 此 这 种 爆发 机 制 称 为 延迟 暴 吕 。 
由 上 述 分 析 ， 我 们 可 以 将 核心 替 缩 型 超新星 爆发 的 整个 过 程 分 为 核心 择 缩 、 产 生 反 弹 激 
= 波 、 激 波 向 外 传播 ， 以 及 原 中 子 星 向 外 辐射 中 微 子 产生 延迟 暴 四 个 阶段 。 在 这 四 个 阶段 中 ， 
P IAA RT Fe 核 的 旋转 、 反 弹 过 程 中 产生 的 对 流 不 稳定 性 和 原 中 子 星 的 振荡 等 都 会 有 相应 的 引 
力 辐射 产生 。 下 面 我 们 将 一 一 介绍 不 同 的 过 程 所 产生 的 引力 波 信 号 。 
2.1 ”旋转 的 Fe AEM Re 
Dimmelmeier 等 人 ”研究 了 旋转 的 Fe Bil AIR S SESERU PEDI Z1, BES 
RATING Fe 核 的 角速度 D.。 在 牛顿 极限 下 ， 较 差 转动 角速度 8 满足 如 下 关系 式 ” : 
A? 
ER PUE. , (1) 
RF, d=rsind 是 到 旋转 轴 的 距离 (0 IRA) SRA 刻画 了 较 差 自转 的 程度 。 当 A 一 oc 
IN, R/R —1, "4 A2 0H, 2/2, = A?/d?. 


表 1 旋转 的 Fe EB ern ^k 93171 S S 0d e 
类 别 ”旋转 速率 2。 最 大 振幅 |hpeax| RAMZ focak HUN HEEL ow 


/(rad- s!) /107” /Hz /1075 Mo- ©) 
1 “fetes c 0.5 700 ~ 800 EIE 
2 6 ~ 13 10 400 ~ 800 5 
3 26~ 13 3.5 e 7.5 70 ~ 200 0.07 ~ 0.5 


Dimmelmeier 等 人 ”计算 了 角速度 2.、 前 身 星 的 质量 M、 参 数 4 以 及 不 同 的 物 态 对 
引力 波 暴 信号 的 影响 。 他 们 的 结果 表明 引力 波 信号 主要 依赖 于 星 核 角 速度 2.， 其 次 才 是 前 
身 星 的 质量 M。 此 外 ， 他 们 发 现 较 差 自转 和 不 同 的 物 态 对 引力 波 的 影响 很 小 ， 具 体 的 计算 


?关于 核心 地 缩 型 超新星 爆发 过 程 的 详细 讨论 ， 可 参看 徐 仁 新 著 的 《天 体 物理 导论 》，2006 年 ， 北 大 出 版 社 ， 第 5.3.3 小 节 。 
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结果 可 参考 文献 Ba 中 的 图 3 和 图 11。Dimmelmeier 及 其 合作 者 根据 Fe 核 角速度 R. 的 大 
小 ， 将 引力 波 信 号 分 为 以 下 三 类 。 

(1) 当 Fe 核 的 旋转 很 慢 时 (BU 2. < 1 ~ 1.5rad. s-!)， 产 生 的 引力 波 的 振幅 将 很 小 ， 其 
最 大 值 [Pasa] (KF 5 x 10 22。 数值 结果 表明 ， 当 缓慢 旋转 时 ， 在 激 波 后 区 域内 将 出 现 强烈 
快速 的 对 流 翻转 ， 这 是 反弹 激 波 速度 降低 后 ， 在 激 波 后 面 的 区 域 形 成 负 箭 梯度 所 致 。 有 具体 
关于 对 流 不 稳定 性 的 问题 我 们 将 在 下 一 节 中 讨论 。 

(2) 当 Q. $6 ~ 13rad- s Ff, Fe 核 旋转 的 速度 加 快 产生 较 大 的 质量 四 极 矩 并 且 内 核 
旋转 的 质量 也 会 增加 ， 因 此 产生 较 强 的 引力 波 信号。 
(3) 24 Fe 核 旋转 得 非常 快 时 ( 即 R. = 13rad . s-!)， 离 心力 作用 的 增加 使 得 核心 反弹 ; 
计算 结果 表明 此 过 程 产生 的 最 大 引力 波 振幅 可 以 达到 [has] = 7.5 x 1077, 

Ott 总 结 了 不 同 的 Fe 核 的 旋转 速度 产生 的 引力 波 信 号 。 在 表 四 中 我 们 展示 了 在 不 
同 的 Fe 核 旋转 速度 下 ， 引 力 波 信号 的 相关 参数 ， 其 中 [hma 是 产生 最 大 引力 波 特 征 振幅 ， 


Ecw 是 辐射 的 总 能 量 ， 广 sx 是 引力 波 能 谱 dBew /df 的 最 大 值 所 对 应 的 频率 。 根 据 Heger 
= 等 人 中 四 关于 恒星 演化 模型 和 脉冲 星 自转 周期 的 计算 结果 ， 有 太阳 金属 丰 度 的 普通 大 质量 
Ea 恒星 的 星 核 旋 转角 速度 低 于 1rad . s-1。 基 于 表 四 的 结果 ， 在 这 种 情况 下 引力 波 的 特征 振幅 
(O 是 0.5 x 1071, SHANA S| AURA BEER ZA 0.1 x 1079 Mo . 2。 若 此 类 源 位 于 银河 系 内 ， 
co 则 该 过 程 产生 的 引力 波 信号 将 会 被 Advanced LIGO 探测 器 探测 到 。 更 快 的 Fe 核 旋 转 的 原 
e 因 可 能 与 伽 马 射线 暴 的 起 源 有 关 ， 相 关 的 研究 可 参看 文献 BI. 

To’ 2.2 ”对 流 不 稳定 性 

© 天 体 物理 所 涉及 的 大 多 数 流体 动力 学 问题 都 会 有 不 稳定 性 。 不 稳定 是 流体 动力 平衡 受 
E 到 某 种 扰动 后 出 现 的 一 种 现象 。 不 稳定 导致 激 波 经 过 后 流体 的 密度 、 压 强 梯度 和 焙 等 的 变 
e 化 ， 促 使 流体 发 生 翻 腾 ， 最 终 使 流体 充分 混合 ， 达 到 平衡 态 温度 ， 这 种 运动 称 为 对 流 不 稳定 
N WE. FERIENE OEN A EB BPRS E PE TE | 79 3 8. 
> 随 着 反弹 演化 的 进行 ， 原 中 子 星 向 外 辐射 中 微 子 ， 将 能 量 传递 给 在 激 波 后 的 净 能 量 沉 


PDK ta, UGA RES E PI PE ABE SH LE BI AE TTL BS BE HES Woosley 和 Weaver" 考 
ET PE IRA ALK BA PR a E EU" 7" (standing-accretion-shock instability, SASI) 
- 产生 的 引力 波 信和 号。 

© Woosley 和 Weaver 研究 了 二 维 的 核心 十 缩 型 超新星 爆发 ， 考 虑 的 前 身 星 的 质量 是 
15 Mo。 整个 模拟 过 程 持续 了 850 ms， 产 生 的 引力 波 特征 振幅 ”可 以 表示 为 : 


1 /2 Gd 
henar( f) — D T? e IJ , (2) 


其 中 ，D 是 波源 的 距离 ，G 是 万 有 引力 常数 ，c 是 光速 。 由 于 中 微 子 驱动 和 长 期 吸 积 激 波 不 
稳定 性 ， 辐 射 的 引力 波 总 能 量 大 约 是 7.5 x 10-12 Ma . cz， 相应 的 频率 范围 是 100 ~ 500 Hz. 
除了 中 微 子 驱 动产 生 对 流 不 稳定 性 以 外 ， 由 于 产生 了 负 的 径 向 轻 子 梯度 ， 原 中 子 星 也 是 对 流 
不 稳定 的 。Miiller 等 人 研究 了 原 中 子 星 的 对 流 过 程 。 在 他 们 的 计算 中 ， 核 心地 缩 型 超 新 
星 爆 发 的 全 过 程 持续 了 1.2 s， 产 生 的 引力 波 的 特征 振幅 是 (2 ~ 5) x 1077, ff S| ZU 


AB 
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的 总 能 量 约 为 1.6 x 1071? Mo .ce2， 特 征 频 率 范围 是 700 ~ 1500 Hz。 我 们 将 相关 结果 展示 在 
口 


表 2 ”对 流 不 稳定 性 过 程 产生 的 引力 波 信号 的 相关 参数 


过 程 特征 振幅 |h| ”特征 频率 f 持续 时 间 At 辐射 能 量 Ew 

/10-?3 /Hz /ms /(107 ? Mo- c?) 
中 微 子 驱动 和 长 期 吸 积 激 波 不 稳定 性 6 100 ~ 500 850 7.5 
原 中 子 星 对 流 2-5 700 ~ 1500 1200 1.6 


3 原 中 子 星 的 振荡 


核心 南 缩 型 超新星 爆发 的 最 终 产 物 是 原 中 子 星 ， 其 内 部 的 温度 很 高 ， 随 着 内 部 的 中 微 子 
向 外 辐射 使 原 中 子 星 内 部 的 温度 逐渐 降低 ， 最 终 原 中 子 星 冷却 。 这 一 章 中 ， 我 们 将 介绍 在 牛 
顿 框架 下 原 中 子 星 的 振荡 理论 ”， 并 分 析 产 生 各 种 振荡 模式 的 恢复 力 ， 然 后 讨论 具体 的 g 模 
式 主导 下 产生 的 引力 辐射 ,最 后 探讨 引力 波 暴 对 中 子 星 物 态 的 依赖 。 

3.1 ”中 子 星 的 非 径 向 振荡 

1988 年 ，McDermott 等 人 "研究 了 在 牛顿 框架 下 中 子 星 的 非 径 向 振荡 。 车 中 子 星 内 部 
某 一 流体 元 发 生 一 小 扰动 后 所 处 的 状态 称 为 扰动 态 ， 此 时 中 子 星 内 部 物理 量 可 以 看 作 是 由 
平衡 态 所 对 应 的 物理 量 加 上 一 小 扰动 变量 构成 。 例 如 ， 密 度 、 压 强 和 引力 势 为 : 


p = potdop 
D = potóp , (3) 
B = +66 


Kh, fo 表示 平衡 态 量 ，6 为 扰动 变量 (fe {p,p,$})。 通 常 描述 流 体 运动 有 两 种 方法 ， 即 
拉 格 朗 日 方法 和 欧 拉 方 法 ， 且 两 种 方法 之 间 存 在 如 下 转换 关系 : 


6f=f'+or-Vfo , (4) 


RP, Of 和 fl 分 别 表示 拉 格 朗 日 扰动 变量 和 欧 拉 扰动 变量 。 
描述 中 子 星 在 线性 、 绝 热 条 件 下 非 径 向 振荡 基本 方程 组 之 前 ，McDermott 等 人 ” 先 将 
MERE Or 分 解 为 径 向 分 量 和 水 平分 量 ， 即 Or = é etén KER, ERRER T, P 
子 星 在 线性 、 绝 热 条 件 下 的 非 征 向 振荡 基本 方程 组 为 : 


连续 性 方程 : | (E) FVn:£j —0 , (5) 
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EREE oE, Of p Oo 100/ 
ETE: 0? pr põr põr’ (6) 
ð én / p 
ETE x VpB Vn ( p ) ) (7) 
eee 0 0 
绝热 状态 方程 : Seu m (8) 
p p 
: . 19 o ; 
泊 松 方程 : r? + Vi) P =4nGo'. (9) 
r? Or Or 
这 里 的 绝热 指数 定义 为 九 三 (djnp/dlnp)aa， 此 外 有 ， 
18 1 ð 1 0? 
= | ; E | 1 
MESSI r 00°? rsind Op Vi ae sin 0 00 (nts : x) r? sin? 0 Oy? (0) 
McDermott 等 人 一 对 此 方程 组 进行 改写 ， 得 到 只 含有 变量 p (7), P(r) 和 总 的 方程 组 : 
1d g S? 1, L(L-4- 1) ， 
sid» (rE) d fr (& 1 dip? PA Pp , (11) 
dp' J 2 2 dg’ 
p — cx ; 12 
dr ^ a p(w -N)6 Par (12) 
1d f ,dp L(L-4- 1) ,, p pN? 
- pe 1 
r? dr (r dr ) r? y ds. (5 H g Sr 13) 


其 中 ，g 三 -~VG， 绝 热 声速 c 满足 c2 = dp/dp = Dp/p FEME (Lamb 频率 ) S; 满足 


S? =1(1+1)c2/r?, AMR (Brunt-Vàisalà 频率 ) N 满足 N? = g ( 


此 外 ， 这 里 还 使 JJ VY (0, p) = m Dy 0p), 其 其 中 Yim(0, 


2 


1 dlap dlnp 
I1 dr dr J’ 


p) 是 球 谐 函数 ，! 是 角 


向 方向 的 节点 数 。 为 了 方便 起 见 ， 我 们 省 略 了 平衡 态 物理 量 的 下 标 “0”。 


A (ILI) — (3) 中 ， 当 径 向 节点 数 [n | 和 水 平方 向 的 节点 数 1 的 值 很 大 时 ， 可 以 忽略 引力 


势 的 扰动 项 ， 即 8' = 0， 这 种 近似 方法 称 为 Cowling 近似 ”。 最 终 式 (中) 一 (03) 改写 为 只 


含有 变量 & 的 二 阶 微分 方程 (详细 的 分 析 过 程 可 参考 文献 EZ, HJ), B 
d£, w? N?\ (8 
dr2 d E s) e: J s Ps 
此 时 ， 式 (T2) 能 够 很 好 地 描述 振荡 模式 的 整体 特性 ， 并 且 也 可 以 得 到 各 种 振荡 模式 所 对 应 
的 频率 范围 。 为 了 能 够 更 好 地 分 析 不 同 的 振荡 模式 ， 我 们 把 式 (a) SH: 
d£, 
as ZTE) > (15) 


Kr) = 5 (5 i) (8 i) (16) 


A. 44= 'HBTII 
11 F RH | | 
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juni 


HX (I3) 我 们 可 以 看 出 ，é&; 的 局 域 表现 形式 依赖 于 天 。(r) 的 符号 。 当 Ks(7) XT 0B. 6. 
可 以 近似 地 表示 为 : 


Er œ cos (frae) j K.(r)>0 ， (17) 


其 中 几 是 由 边界 条 件 所 决定 的 相位 。 式 (I2) AH ERA, 6, 是 7 的 振荡 函数 。 当 
K,(r) NF ORT, £ 可 以 近似 地 表示 为 : 


En & exp (= / IK, [ar] K,(r) <0. (18) 

由 式 (IN) 我 们 可 以 看 出 ，é&; 是 7 的 指数 递增 或 递减 的 函数 。 这 样 ， 我 们 根据 w 的 取 值 分 为 
两 种 情况 ， 

lu] >|N| 和 Jw|»5 , (19) 

lwl <N] 和 jul «s . (20) 


相应 地 ， 我 们 可 以 得 到 两 类 振荡 模式 : 满足 式 (四 ) 的 振荡 模式 称 为 高 频率 振荡 模式 ， 其 恢 
复 力 是 压力 ， 所 以 高 频率 振荡 模式 又 称 为 p 模式 ; 满足 式 (ZU) 的 振荡 模式 称 为 低频 率 振荡 
模式 ， 其 恢复 力 是 浮力 ， 因 此 将 低频 率 振荡 模式 称 为 9 模式 。 此 外 ， 与 p 模式 类 似 的 一 种 振 
荡 模 式 是 上 模式 。 虽 然 f 模式 的 恢复 力 仍 然 是 压力 ， 但 其 频率 范围 却 介 于 p 模式 和 9 模式 
之 间 。 关 于 f 模式 的 详细 讨论 可 参考 McDermott 5 A" 的 研究 工作 。 更 多 的 恒星 振荡 模式 
可 参考 文献 A, bU]. 
3.2” 原 中 子 星 g 模式 主导 下 产生 的 引力 波 暴 

McDermott 等 人 ”计算 了 在 线性 、 绝 热 条 件 下 ， 原 中 子 星 的 各 种 振荡 模式 产生 的 
引力 辐射 以 及 各 种 振荡 模式 所 对 应 的 振荡 周期 。2006 4E, Ott 等 人 "研究 了 二 维 核心 地 缩 
型 超新星 在 核心 南 缩 、 反 弹 阶 段 以 及 原 中 子 星 的 9 模式 下 产生 的 引力 波 暴 。 他 们 选择 了 
s11WW"', m15b6"" f s25WW 三 种 模型 ， 其 前 身 星 的 质量 分 别 为 11 Mo, 25 Mo, 15 Mo: 
此 外 ，s11WAW 模型 、m15b6 模型 和 s25 WW 模型 中 Fe 核心 质量 分 别 是 1.37 Mo, 1.47 Mo 
和 1.92 Ms. Ott 等 人 号 通过 计算 发 现 ， 对 于 sl. WW 模型 和 m15b6 模型 ， 当 演化 时 间 持 续 
到 400 ms 时 ， 引 力 波 完全 由 9 模式 主导 产生 。 通 过 Ott 等 人 ™ 的 计算 结果 ， 我 们 发 现在 整 
~ 个 超新星 爆发 的 过 程 中 ，g 模式 主导 产生 的 引力 波 是 最 强 的 ， 与 之 相对 应 的 引力 波 特 征 振幅 
也 最 大 。 而 对 于 s25WW 模型 ， 其 前 身 星 的 星体 结构 与 前 两 种 模型 相 比 有 很 大 的 差异 ， 更 
重要 的 是 s25WW 模型 中 Fe 核心 的 质量 大 于 前 两 种 模型 ， 因 此 在 演化 的 过 程 中 产生 的 引力 
波 波形 也 与 前 两 种 模型 不 一 样 。 在 s25WW 模型 中 ，g 模式 主导 产生 的 引力 波 发 生 在 500 ms 
时 (具体 的 引力 波 波形 图 可 参看 Ott 5 AP" 研究 工作 中 的 图 1 和 图 2)， 当 演化 时 间 持 续 到 
900ms 时 ， 产 生 了 更 强 的 引力 波 信号 。Ott 等 人 分 析 认 为 ， 此 时 的 引力 波 信 号 是 由 更 高 
阶 的 9 模式 产生 ， 其 引力 波 特 征 振幅 可 以 达到 5 x 10-?0。 

表 回 展示 了 原 中 子 星 9 模式 产生 的 引力 波 信号 的 相关 参数 ， 其 中 特征 振幅 Rehar 由 公式 
加 给 出 。 从 表 田 中 我 们 可 以 看 到 ，s11WW 模型 和 m15b6 模型 产生 的 引力 波 特征 振幅 约 为 


ChinaXiv 合 作 期 刊 


518 天 文学 进展 39 卷 


1.3 x 10?!, m S25 WW 模型 产生 的 引力 波 特征 振幅 可 以 达到 5 x 10-20。 此 外 ， 三 种 模型 中 
辐射 的 引力 波 的 总 能 量 分 别 是 1.4 x 1078 Mo -c?, 1.6 x 1078 Mo -2 7 8 x 1075 Mo c?. He 
据 表 回 展 示 的 引力 波 信号 的 相关 参数 ， 我 们 发 现 其 前 身 星 的 质量 和 Fe 核心 质量 越 大 ， 产 生 的 
引力 波 的 能 量 也 越 大 。Ott” 通过 不 同 的 模型 计算 发 现 ， 在 整个 核心 韦 缩 型 超新星 的 过 程 中 
中 子 星 的 振荡 模式 可 能 是 产生 引力 波 的 主要 机 制 |。 


X3 g 模式 主导 下 产生 的 引力 波 信号 的 相关 参数 加 


模型 特征 振幅 |h| ”特征 频率 f ”持续 时 间 At 辐射 总 能 量 Ew 
/107! /Hz /ms /(10-7 Mo :©2) 
s11WW 1.3 654 1045 0.16 
s25WW 50.0 937 1110 824.28 
m15b6 1.2 660 927.2 0.14 


3.3 ”引力 波 暴 的 物 态 依赖 性 

最 后 我 们 讨论 引力 波 暴 的 物 态 依赖 性 。 我 们 知道 ， 物 态 方程 是 描述 星体 微观 属性 和 宏观 
结构 的 重要 纽带 。 物 态 方程 给 出 了 物质 内 部 压强 p 和 密度 p 之 间 的 关系 ， 如 果 密 度 增 加 所 对 
应 的 压强 上 升 较 快 ， 则 称 物 态 方程 较 “ 硬 ”; 反之 ， 则 较 “ 软 ” 研究 结果 表明 ， 不 同 的 星体 
内 部 成 分 和 结构 ， 甚 至 同一 种 结构 中 运用 不 同 的 相互 作用 理论 ， 都 会 给 出 不 同 的 物 态 方程 。 

Marek 等 人 名利 用 二 维 情况 下 质量 为 15 Mo 的 前 身 星 模 型 人 ~"， 模 拟 了 核心 地 缩 型 超 新 
星 爆 发 的 全 过 程 。Marek 等 人 在 研究 过 程 中 采用 了 两 种 不 同 的 核 物 态 方程 ， 即 Lattimer 和 
Swesty” 物 态 方程 (L-S 模型 )， 以 及 Hillebrandt 等 人 ”的 物 态 方程 (H-W 模 型 )。 两 种 模型 
最 大 不 同 就 是 H-W 物 态 方 程 比 L-S 模型 更 硬 。Marek 等 人 计算 发 现 ，L-S 模型 产生 的 引力 
波 特 征 振 幅 是 5 x 10-22， 特 征 频率 范围 在 600 ~ 800 Hz 内 。 与 L-S 模型 相 比 ，H-W 模型 
计算 得 到 的 引力 波 特征 振幅 更 小 ， 相 应 的 特征 频率 范围 在 300 ~ 600 Hz 内 。 计 算 结 果 表 明 ， 
L-S 模型 物 态 更 “ 软 ”， 中 子 星 致密 度 更 大 ， 产 生 的 引力 辐射 也 更 强 。 


P i 


c 4 总 结 与 展望 


我 们 在 本 文 详细 讨论 了 Fe 核 旋转 的 角速度 和 对 流 不 稳定 性 对 引力 波 暴 的 影响 ， 并 且 分 
析 了 中 子 星 的 各 种 振荡 模式 ， 包 括 g, f, p 模式 。 在 此 基础 上 ， 我 们 讨论 了 9 模式 产生 的 引 
力 波 信号 的 特征 振幅 、 频 率 以 及 辐射 的 总 能 量 。 与 Fe 核 旋 转 的 角速度 和 对 流 不 稳定 性 相 比 ， 
引力 波 主要 来 源 于 原 中 子 星 的 振荡 模式 。 除 上 述 振荡 模式 外 ，Andersson 和 Kokkotas"" fF 
究 发 现 ， 若 不 考虑 中 子 星 内 部 的 耗 散 机 制 ， 转 动 的 中 子 星 的 7 振荡 模式 是 不 稳定 的 。 在 此 基 
fili E, Owen 和 Lindblom"" 研究 发 现 ，" 模式 的 不 稳定 性 产生 的 引力 波 信号 将 有 可 能 被 新 一 
代 的 引力 波 探测 器 探测 到 。 

虽然 核心 替 缩 型 超新星 是 重要 的 引力 波 天 体 物 理 源 之 一 ， 但 是 至 今 地 面 引 力 波 探测 器 都 
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没有 直接 探测 到 相关 的 引力 波 信 号 ””， 其 原因 是 超新星 爆发 时 复杂 的 物理 过 程 以 及 中 子 
星 的 物 态 至 今 尚 未 清楚 。 因 此 ， 即 使 理论 上 有 相应 的 波形 模板 ， 也 未 必 精 确 。 而 GW170817 
引力 波 事件 ”提供 了 解决 这 些 问 题 的 办 法 。 该 引力 波 事件 是 由 相互 绕 转 的 两 个 中 子 星 旋 近 阶 
段 产 生 ， 更 重要 的 是 有 多 波段 的 电磁 辐射 观测 。 因 此 对 于 核心 去 缩 型 超新星 爆发 ， 我 们 亦 可 
采用 相同 的 方法 ， 对 这 类 事件 进行 电磁 波 和 引力 波 的 协同 观测 ， 这 样 将 会 更 好 地 理解 核心 坪 
缩 型 超新星 爆发 的 整个 过 程 ， 从 而 进一步 来 探究 引力 波 暴 的 波形 特征 。 

最 后 值得 一 提 的 是 ，2020 Æ 12 H 10 日 4 时 14 分 ， 由 中 国 科 学 院 空间 科学 (二 期 ) 
先导 专项 部 署 的 引力 波 暴 高 能 电磁 对 应 体 全 天 监测 器 (gravitational wave high-energy 
electromagnetic counterpart all-sky monitor, GECAM) 卫星 顺利 发 射 升 空 ”。GECAM E 
星 的 科学 目标 是 对 引力 波 伽 马 暴 、 快 速射 电 暴 高 能 辐射 等 高 能 天 体现 象 进行 全 天 候 的 观测 。 
我 们 也 和 希望 未 来 通过 GECAM 与 地 面 引 力 波 探测 器 的 协同 观测 ， 能 够 深入 研究 核心 替 缩 型 
超新星 的 爆发 过 程 以 及 中 子 星 的 内 部 结构 。 
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Abstract: Ever since the very first experimental effort to detect gravitational waves, core- 
collapse supernovae have been considered as important gravitational wave sources. In this 
paper, we review the relevant knowledge of gravitational wave in general relativity. We 
introduce the influence of the angular velocity of the iron core and the convection instability 
on gravitational wave bursts in stellar core collapse and its post-bounce phase. We overview 
the oscillation theory of proto-neutron stars, the final product of core-collapse supernova 
explosion, and the resulting gravitational wave bursts. Finally, we expect that the coincident 
observation of gravitational waves and electromagnetic waves by the international network 
of observatories, can help us understand the core-collapse supernova evolution and internal 


structure of neutron stars in the near future. 
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